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ABSTRACT 



The quantization of a constant of motion for the harmonic oscillator with a time- 
explicitly depending external force is carried out. This quantization approach is com- 
pared with the normal Hamiltonian quantization approach. Numerical results show that 
there are qualitative and quantitative differences for both approaches, suggesting that 
the quantization of this constant of motion may be verified experimentally. 
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I. Introduction 



Modern Physics has emerged as theories formulated in terms of Hamiltonians or 
Lagrangians structures. Among them, one could mention classical physics, field theories 
[1], quantum mechanics [2] and statistical mechanics [3]. Nature has been kind enough 
to allows its description , so far, in terms of Hamiltonian or Lagrangian formulation, 
despite of some mathematical problems that these formulations have by themselves [4]. 
These problems appear mainly for dissipative [5] or time-explicitly depending phenomena 
[6], where one wonder whether the associated constant of motion or Hamiltonian is the 
relevant parameter to study the associated physics. For the harmonic oscillator with a 
time-explicitly depending force (HOTED), it is possible to know a constant of motion 
which is also time-explicitly depending [7]. Therefore, one may wonder about any differ- 
ence that HOTED could have from the quantum mechanics point of view. In this paper, 
the quantization of a constant of motion fort HOTED is associated and compared with 
the usual Hamiltonian approach. 



II. Constant of Motion and Quantization Approach 

Consider that the harmonic oscillator is perturbed by the time-explicitly depending 
force given by 

f(t) = Asmttt, (1) 

where A and O are the amplitude and the frequency of the force. The resulting classical 
dynamical system can be written as 

dx . . 

*=" <2a) 



and 



— = — u x H sin Sit , (26) 

at m 



where u is the frequency of the free harmonic oscillations. A constant of motion associated 
to [2] for O / u is given by [7] the function 



m A ( \ A 2 

K(x,v,t) = — (v 2 +u 2 x 2 ) + ^ o [ ftv cos fit +u 2 x sin tit * sin 2 (Ot) . 

2 ft 2 - oj 2 V / 2m(0 2 - u 2 ) 

(3) 



One can readily see that (3) is a constant of motion since it satisfies the following partial 
differential equation 

dK o A . ^ , dK dK 

v — + [-u 2 x + — sin fit — + — = . 4 

ox m ov Ot 

The approach used here to quantize (3) is based on the construction of the associated 
Shrodinger's equation 

ih—=K{x,v,t)V , (5) 

where \I/ is the wave function, ^ = \I/(;r,t), and K is the hermitian operator associ- 
ated to function (3). This hermitian operator is constructed by making the following 
substitutions 

x — >■ x = x and v — >■ v = — i — — (6) 

m ox 

which bring about the Shrodinger's equation 

fl* Ti 2 dH 1 2 2 T 

in-p— = — 7T + -mw x W 

dt 2m dx 2 2 

A 2 2 ^ Atu 2 x 

or sin fit H o ■ 

2m(fi 2 - u 2 ) fi 2 - u' 

. Am cos fit d^ 
m(fi 2 - a; 2 ) dx 

The first line on Eq. (7) represents the pure harmonic quantum oscillator. The eigen- 
functions of the pure harmonic oscillator are well known [8] and are better handle with 
annihilation and creation operators, a and a + . x and v are written in terms of these 
operators as 

(a + a + ) , (8a) 



+ \ ~7> — ~ 27 sm 2 0t+ 7 ^ 2 sin fit \ * (7) 



x = 



2mw 



and 



J^(<> <>')■ (86) 



Then, Eq. (3) can be written as 

yl 2 

— ^- sin^ S, 2t 
- u z \ 

(9a) 



X = K (a + , a) - n 2 27 sin2 n * 

2m(il — to ) 



H o ?\/— {[a; sin fit — ificosfitla + [a; sin fit + ificosfit|a + 

fi 2 - a; 2 V 2m I ' 



where K (a + 1 a) is given by 

K = hu{a + a + 1/2) . (96) 

If |n > represents an eigenfunction of the pure harmonic oscillator, K a , one has the usual 
properties 

(10a) 



a\n > \pn \n — 1 > , 
a + \n >= \/n + 1 |n + 1 > 
a + a|n >= nln > 



and 



(106) 
(10c) 

(10c) 



[a,a + ] = 1 , 

where symbol [, ] represents the commutator between operators, [a,a + ] = aa + — a + a. 
Now, proposing in (5) a solution of the form 



V(x,t) = J2 c n(t)\n > 



(11) 



n=0 



one gets the system of equations 

A 2 



ihc r 



+ 



hu(n + 1/2) - 
A 



2m(ft 2 - u 2 ) 



sin 2 fit 



I hu^ | ^ — — - ^ ^ ^ ^ Ot)c m+1 + \fm{uj sin fit + ifi cos fit) c r 
V 2m [ 



n 2 -u 2 

which can be written in terms of the parameters 

r = cut , e = A/hfl, fi = Ti/mtjj and p = fi/w 

as 



(12) 
(13) 



zcL = ( n + 1/2 - 27 sin 2 pr I c m 



2ft(l - p 2 ) 

+ e\/2h (1 — p) | \Jm + l(sin pr — zp cos pr)c m+ i + \/™(sin pr + z'p cos pr)c m _i j , 

(14) 

where c^ denotes the expression d m = dc m /dT. 

On the other hand, the Hamiltonian quantization approach for the system (2) is 
formulated by the Shrodinger's equation 



where the Hamiltonian operator H is given by 

•~ v 2 1 

H = — + -mu 2 x 2 - xA sin fit . (16) 
2m 2 v 1 

The operators £ and p can be written in terms of the annihilation and creation operators, 
a and a + , as 

i = V&° + a+) (17o) 

and 



0=-^ (.-„+), (176) 

where a and a + satisfy the same relations (10). Written (16) in terms of these operators, 
it follows that 

H = H a (a + ,a) - J— — Asinttt (a + a + ) , (18) 
V 2mu) 

where H a (a + ,a) is given by H Q = Tiu(a + a + 1/2). Using the same expansion (10) and 
parameters (13), one gets the system of equations 

ic' m = (m + l/2)c m - epsinpr {y/m c m -i + Vm+ 1 c m +i) ■ (19) 



III. Numerical Solution 

To know the evolution of c m (r) for Eqs. (14) and (15), these can be solve numerically. 
The fixed parameters used are ft = 0.4 and p = 6.25, and one studies the dependence of 
the system on the parameters e and r. The initial conditions are 

c o (0) = l, ,c,(0) = for J = 1,2,... (20) 

The Fig. 1 shows the evolution of the probabilities |c (r)| 2 and |ci(r)| 2 and their 
dependence on the parameter e. The upper plots correspond to the quantization of the 
constant of motion (Eq. 14), and the lower plots correspond to the quantization of the 
Hamiltonian (Eq. 19). For the same e, the peak values of the probabilities |c (r)| 2 and 
|ci(t)| 2 occur at different times (r) and with different amplitudes. 

The Fig. 2 shows the maximum number of exited states involved in the dynamics 
of the system for a given e. The criterion used is the following: one defines that a state 
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\n > is involved in the dynamics if the probability that the system to be in this state is 
higher than 0.0001, \c n (r)\ 2 > 0.0001. As one can see, for a given e the number of levels 
involved in the Hamiltonian quantization can be much larger than the levels involved in 
the constant of motion quantization. 

The Fig. 3 shows the evolution of the expected value of x 2 7 

oo 

< x 2 >=< *|x 2 |* >= <4(T)c n (V) < m\x 2 \n > , (21) 

m,n=0 

for the constant of motion and Hamiltonian quantization. As one can see, the value 
for the Hamiltonian case is about one order of magnitude higher than its value for the 
constant of motion case. 



III. Conclusions 

The quantization of the harmonic oscillator with a time-explicitly depending external 
force has been studied from the point of view of the constant of motion. The results have 
been compared with the known case of Hamiltonian approach. The probabilities |c G (r)| 2 
and |ci(t)| 2 , the expected value < x 2 > and the number of exited states involve in the 
dynamics are quite different for the constant of motion and Hamiltonian approaches. 
These results may be checked experimentally taking care that the external force added 
to the harmonic oscillator system is of the form given by Eq. (1). 
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Figure Captions 

Fig.l Evolution of the probabilities |c (r)| 2 and |ci(r)| 2 for 

the system (14) and (19). The value of the other parameters are 

h = 0.4 and p = 6.25. The initial conditions are given by Eq. (20). 

Fig. 2 Number of exited states involved in the dynamics as a function of the 
parameter e. The other parameters are the same as Fig. 1. 

Fig. 3 Evolution of the expected value < x 2 > for the constant of motion and 
Hamiltonian quantization approaches. The parameter epsilon takes the values 
e = 0, 5, 10. The other parameters and initial conditions are the same as 
Fig. 1. 
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